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ABSTRACT 

We present numerical simulations of the evolution of low-mass, isothermal, molec- 
ular cores which are subjected to an increase in external pressure P cx t- If -P ex t increases 
very slowly, the core approaches instability quite quasistatically. However, for larger 
(but still quite modest) dP ey ^/dt a compression wave is driven into the core, thereby 
triggering collapse from the outside in. If collapse of a core is induced by increasing 
P e xti this has a number of interesting consequences, (i) The density profile is approxi- 
mately flat in the centre during the prestellar phase (i.e. before the compression wave 
converges on the centre creating the central protostar) . (ii) During the prestellar phase 
there are (subsonic) inward velocities in the outer layers of the core, whilst the inner 
parts are still approximately at rest, (iii) There is an initial short phase of rapid accre- 
tion (notionally the Class phase) , followed by a longer phase of slower accretion (the 
Class I phase). All these features accord well with observation, but are at variance 
with the predictions of the standard theory of star formation based on the inside-out 
collapse of a singular isothermal sphere. We note that the setting up of a coherent 
inward velocity field appears to be a generic feature of compression waves; and we 
speculate that interactions and interference between such velocity fields may play a 
crucial role in initiating the fragmentation of cores and the genesis of multiple star 
systems. 
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1 INTRODUCTION 

The initial conditions for star formation are in general rather 
weakly constrained. The observations (e.g. Andre, Ward- 
Thompson, & Barsony 2000; Myers, Evans, & Ohashi 2000) 
are limited by telescope resolution and confusion, and their 
interpretation is hindered by the complexities and uncer- 
tainties of gas-phase abundances, excitation conditions, and 
radiation transport. However, a picture is emerging which 
suggests that star formation is sometimes triggered rather 
impulsively. 

Molecular-line mapping by Myers and collaborators 
(e.g. Myers & Benson 1983; Benson & Myers 1989) and sub- 
sequent searches by Clemens & Barvainis (1988), Bourkc, 
Hyland & Robinson (1995), Jessop & Ward- Thompson 
(2000), have established a large sample of dense cores which 
appear to be the sites of ongoing or imminent star formation. 
Many of these cores contain IRAS sources (Beichman et al. 
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1986), and are therefore presumed to have already formed 
protostars. Those which do not contain IRAS sources are 
termed starless cores. 

Submillimeter mapping of these starless cores by Ward- 
Thompson et al. (1994, 1999) has identified a subset which, 
on the basis of their virial ratios and relatively high cen- 
tral densities, are likely to be in a state of imminent or on- 
going contraction. These starless cores are therefore termed 
prestellar (strictly pre-protostellar) . Submillimeter mapping 
of prestellar cores (Ward-Thompson et al. 1994, Andre et 
al. 1996, Ward- Thompson et al. 1999) suggests that, if the 
emitting dust is isothermal, the cores have rather flat den- 
sity profiles in the centre. Specifically, r\ = —d£n[p]/d£n[r] is 
in the range to 1, in the innermost few thousand AU. This 
conclusion may be somewhat weaker if, as seems likely (Jes- 
sop & Ward-Thompson 2001; Evans et al. 2001; Zucconi, 
Walmsley & Galli 2001; Ward- Thompson, Andre & Kirk, 
2001), the dust temperature decreases towards the centre of 
a core. However, mid-infrared absorption measurements of 
prestellar cores - which are insensitive to the temperature 
profile - also indicate flat inner density profiles (Bacmann 
et al. 2000) . Additionally, mid-infrared observations suggest 



2 P. Hennebelle, A. P. Whitworth, P. P. Gladwin & Ph. Andre 



that beyond ~ 10,000 AU the density profile may steepen 
to t) ^ 4 (Abergel et al. 1996, Bacmann et al. 2000). This 
steepening of the outer envelope appears to occur at smaller 
radii in relatively close-packed protoclusters like p Ophiuchi 
than in regions of distributed star formation such as Taurus; 
viz. at ~ 3,000 AU in p Ophiuchi, and at ~ 15,000 AU in 
Taurus (Motte & Andre 2001). 

The asymmetric self-absorbed molecular-line profiles 
observed in some prestellar cores (Tafalla et al. 1998, 
Williams et al. 1999, Lee, Myers & Tafalla 1999, Gregersen & 
Evans 2000) suggest that they are indeed already collapsing. 
In particular, the detailed analyses of L1544 by Tafalla et al. 
(1998) and Williams et al. (1999) imply that the inner parts 
of the core are relatively stationary, and an approximately 
uniform velocity field has been established in the outer lay- 
ers. This is very reminiscent of the velocity fields which are 
set up by inward-propagating compression waves in similar- 
ity solutions for contracting isothermal spheres (Whitworth 
& Summers 1985). It is this apparent similarity which we 
explore in the present paper. 

The prestellar phase terminates as soon as a star-like 
object forms at the centre of the dense core, and the core 
then becomes a Class protostar. Conceptually, the Class 
protostellar phase terminates once the extended envelope 
contains less than half the total mass of the original core, i.e. 
more than half the mass is in the central star-like object(s) 
plus attendant disc(s) (Andre, Ward- Thompson & Barsony 
1993, 2000); the Class I phase then begins (Lada & Wilking 
1984, Lada 1987, Andre & Montmerle 1994). The Class I 
protostellar phase terminates when most of the envelope has 
been accreted or dissipated, revealing a classical T Tauri star 
(CTTS) accreting from a residual circumstellar disc, i.e. a 
Class II object. Once the inner disc has been dissipated, the 
accretion rate is greatly reduced and the source becomes a 
weak-lined T Tauri star (WTTS) or Class III object. 

On the basis of statistical arguments (i.e. source num- 
bers and a presumed constant star formation rate) it is in- 
ferred (Beichmann et al. 1986, Andre, Ward- Thompson & 
Barsony 2000) that the prestellar phase lasts, typically, 10 6 
to 10 7 years. In close-packed protoclusters like p Ophiuchi, 
the Class phase appears to last a few times 10 4 years, and 
is characterized by powerful collimated outflows indicative of 
rapid accretion, <; 10~ 5 M Q year -1 (Bontemps et al. 1996). 
In distributed star formation regions like Taurus, the dura- 
tion of the Class phase appears to be longer, ~ 10 5 years, 
and the accretion rates lower, ^2x 10 -6 Mq year" 1 (Motte 
& Andre 2001). The Class I phase appears to last ~2x 10 5 
years (e.g. Greene et al. 1994, Kenyon & Hartmann 1995) 
and is characterized by slower accretion, ^ 10~ 6 Mq year -1 , 
and weaker less collimated outflows (Bontemps et al. 1996). 
Together, the Class II and Class III phases appear to last 
<; 10 7 years, but seemingly the transition from Class II 
(CTTS) to Class III (WTTS) is very short and can occur at 
any time; there are both very young WTTSs (close to the 
birthline, on the right of the Hertzsprung-Russell Diagram) 
and very old CTTSs (approaching the Main Sequence on 
the left of the Hertzsprung-Russell Diagram) (e.g. Stahler & 
Walter 1993). 

These observationally inferred features of prestellar and 
protostellar evolution combine to form a reasonably coher- 
ent picture. However, especially in star-forming clusters, 
this picture is difficult to reconcile with the standard the- 



ory of Shu, Adams & Lizano (1987) based on the inside- 
out collapse of a singular isothermal sphere. The strength 
of the standard theory is that it makes specific quantita- 
tive predictions, but some of these predictions are difficult 
to reconcile with observation, (i) In the standard model, 
prestellar cores should be strongly centrally condensed, r) = 
—d£n[p]/d£n[r] ~ 2 and young protostars should be some- 
what less centrally condensed, rj ~ 3/2 . By contrast, obser- 
vations suggest that prestellar cores have rather flat central 
density profiles, r\ Si 1, and Class protostars have steeper 
ones, (ii) The standard theory predicts that prestellar cores 
are static, and that inward motions only develop during the 
protostellar phase and are initially confined to the central 
regions. In contrast, the observations imply that inward mo- 
tions already exist during the prestellar phase, and that ini- 
tially they are more rapid in the outer regions, (iii) The stan- 
dard theory predicts that the accretion rate is roughly con- 
stant, and hence the Class and Class I lifetimes should be 
comparable. Observations detect many more Class I sources 
than Class ones - although this statistical result may be 
compromised by the difficulty of measuring precisely when 
an object has accreted half the total mass of its initial core, 
(iv) Additionally, the standard theory assumes initial con- 
ditions which are unlikely to arise in nature, because they 
are both singular and unstable, (v) Also the standard theory 
has a strong inbuilt pre-disposition to the formation of single 
stars - in stark contrast with the high proportion of bina- 
ries and higher multiples observed in young star-formation 
regions. Therefore, although the standard theory may pro- 
vide a good, zeroth-order description of protostellar collapse 
in sparse, quiescent star-formation regions like Taurus (cf. 
Motte & Andre 2001), it appears that more dynamical mod- 
els are required to understand close-packed regions like p 
Ophiuchi. 

Foster & Chevalier (1993) have explored how the col- 
lapse of an isothermal core develops if one abandons the 
assumption of singularity. Their simulations start from non- 
singular isothermal equilibria, i.e. Bonnor-Ebert spheres 
(e.g. Bonnor 1956), and collapse is then triggered by discon- 
tinuously increasing the density. The ensuing collapse results 
in supersonic inflow velocities in the centre of the core at the 
end of the prestellar phase, and in a marked decline of the 
accretion rate from the Class to the Class I phase (see 
Henriksen, Andre, Bontemps 1997 and Whitworth & Ward- 
Thompson 2001 for simpler, pressure-free descriptions of this 
evolution) . 

In this paper we pursue the consequences of non- 
singularity further. Our simulations also start from non- 
singular isothermal equilibrium cores, but collapse is then 
triggered by a steady increase in the external pressure (cf. 
Myers & Lazarian 1998). Most of the observational con- 
straints detailed above can be reproduced rather well by 
this model. 

In Section 2 we describe the numerical method we use, 
and the initial and boundary conditions. In Section 3 we 
present the results, and in Section 4 we discuss them. Section 
5 summarizes our main conclusions. 
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2 NUMERICAL METHOD 
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The simulations are performed using a Smoothed Particle 
Hydrodynamics code - essentially that described in Turner 
et al. (1995). We use the B2-spline smoothing kernel of Mon- 
aghan & Lattanzio (1985), 

r 4 - 6s 2 + 3s 3 , < s < 1 ; 
W(s) = —{ (2~sf, Ks<2; (1) 

47r I 0, s>2. 

and the smoothing lengths of the individual particles are ad- 
justed so that the kernel encompasses A/" ne ib — 50 ± 5 neigh- 
bours, i.e. for particle i the smoothing length hi is adjusted 
so that there are A/" ne ib other particles j having 
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then the equation of motion for particle i is 
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The first summation determines the hydrostatic and 
viscous accelerations. Since W'(s) = for s > 2 , this sum- 
mation only extends over neighbours, i.e. particles satisfying 



Eqn. 



Ilij- is the artificial viscosity term and acts only 



between neighbours which are approaching one another. We 
use a = 1 , /3 = 2 , and 7 = 0.01 . 

The second summation determines the gravitational ac- 
celeration. Since W*(s) > for all s > , this summation 
extends over all particles. The gravitational attraction be- 
tween particles whose kernels overlap is softened by W*(s) , 
in accordance with Gauss's Gravitational Theorem; implic- 
itly the gravitational softening length and the hydrostatic 
smoothing length are the same, as advocated by Bate & 
Burkert (1997). We use an octal spatial tesselation tree 
(Barnes & Hut 1986, Hernquist 1987) to find neighbours, 
and also to speed up the calculation of gravity. Thus in the 
second summation in Eqn. (^) j' may represent the identifier 
of an individual particle or a cell in the tree. To evaluate this 
sum, the tree is walked from the root cell (representing the 
whole computational domain) downwards. Whenever a cell 
j' is encountered which satisfies the non-opening condition 



that cell contributes to the sum as a point mass, and all 
smaller cells and individual particles within that cell can 
be neglected. In Eqn. (B), Lji is the linear size of the cell 
under consideration, R^/ is the position of the cell's centre 
of mass, and 9 C ii% is the maximum opening angle which a 
cell can subtend at particle i without being opened. We use 
6> crit = 3" 1/2 ~ 0.577. 

The density at particle i is given by 
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We use an isothermal equation of state 
di — a a , Pi — a a pi , 

and so we do not need to solve an energy equation. 
The code uses multiple-particle time-steps, 
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The maximum possible time-step for each particle is given 
by 
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Then each particle i is allocated the largest At* 
smaller than At™ ax , i.e. At Hi with 



(9) 



which is 



m 



INT 



■[Aff""7At ] 
£og[2] 



(10) 



with the additional proviso that the current time-step must 
be an integer multiple of At ni ■ 

At the start of a simulation, cores are modelled as trun- 
cated equilibrium isothermal spheres, contained by a hot 
rarefied external medium. These initial conditions are real- 
ized by distributing particles randomly in a cubic box with 
periodic boundary conditions, and settling them using SPH 
forces only (i.e. the terms in the first summation in Eqn. (^)); 
this creates a uniform-density settled cube. Next a sphere is 
cut from the settled cube, and the particle positions within 
the sphere are stretched radially to reproduce the density 
profile of a truncated isothermal sphere. Isothermal spheres 
having different degrees of central condensation are created 
by picking different values of £b = Ro (4-7rGp c ) 1//2 a^ 1 , where 
£b is the dimensionless radius of a truncated equilibrium 
isothermal sphere, R is the physical radius of the core (not 
the scale- length), and p c is the central density in the core. 
Stability requires £b < 6.45 (Bonnor 1956). 

The particles inside the sphere are then tagged as be- 
ing internal particles, and henceforth they experience both 
SPH forces and mutual gravity, i.e. they are evolved us- 
ing both the summations in Eqn. (^). The particles outside 
the sphere are tagged as external particles, and henceforth 
they are evolved using only SPH forces, i.e. only the first 
summation in Eqn. (W). The external particles are given 
whatever sound speed is needed for them to deliver the pre- 
scribed external pressure P e xt(t) at the edge of the sphere, 
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i.e. tti(i) = (P cx t(t)/pi(i)) 1 ^ 2 • In addition, there is a layer 
of fixed edge particles at the boundary of the computational 
domain. These fixed edge particles have constant density 
pedgc and they are not evolved, but they exert SPH forces 
on the external particles and their sound speed is given by 

ledge = (-Pext(i)/Pedgc) 1 ^ 2 • 

At the outset, the sphere is relaxed with P ex t held con- 
stant, so that it can settle to equilibrium. There is some 
variance between the true isothermal sphere characterized 
by the chosen value of £ D and the actual SPH-modelled con- 
figuration, due to the density smoothing at the core bound- 
ary, where the density gradient is steep. This leads to an 
effective external pressure which is a little higher than the 
specified one. The only significant consequence of this is that 
the onset of collapse is somewhat accelerated in the case of 
slow compression (<f> — 10). This is due to transient oscil- 
lations which for slow compression have time to propagate 
inwards from the boundary and increase the density at the 
centre. 

Stable equilibrium isothermal spheres constructed in 
this way are then subjected to increasing external pressure. 
For simplicity we prescribe 



Pext(t > 0) = P cxt (0) + Pot, 



(11) 



with constant P . Rather than give P in physical units, we 
define a dimensionless parameter 



Pcxt(0)/Po 
Po/lo 



(12) 



which is the ratio of the time-scale on which the external 
pressure doubles to the initial sound- crossing time (P is the 
initial radius of the core) . Small P corresponds to large <j> . 
Simulations have been performed with a variety of <j> values 
between <j> = 10 (very slow compression) and 4> — 0.1 (very 
fast compression). 

There is a simple spherical sink at the centre of the core, 
having radius 150 AU (i.e. ~ 1% of the initial core radius), 
and representing the central protostar (plus attendant disc). 
Any particle entering the sink is assimilated by it. The pro- 
tostar is assumed to come into existence as soon as the sink 
starts to gain mass. Hence we can easily compute its mass 
and accretion rate, as functions of time. 

We fix the core mass at M — 1 Mq , the sound speed 
in the core at a — 0.2 km s -1 (corresponding to a cosmic 
mixture of molecular hydrogen, helium and heavy elements 
at ~ 10 K), and the stability parameter of the initial core 
at £ D = 3.0. The radius of the core is initially ~ 0.08 pc. 
The standard simulations presented in the next section are 
performed with A/int ~ 50,000 internal particles, A4xt ~ 
25,000 external particles, and A/idge ~ 20,000 edge particles; 
some simulations with fewer particles have been performed 
to demonstrate convergence. 



3 RESULTS 

Figure jj] shows detailed results for the <^> = 10 case (i.e. 
very slow compression; P ext takes ten sound-crossing times 
to double). In this case the evolution towards gravitational 
instability is relatively quasistatic, because there is time 
for sound waves to slop around re-distributing the mat- 
ter. During the prestellar phase (the first three rows on 



Fig. |lj), the outer boundary is pushed inwards, and the 
core becomes increasingly centrally condensed; the outer 
parts of the core move slowly inwards at speeds up to 
~ 0.03kms _1 (= 0.15a o ), whilst matter near the centre is 
approximately stationary. Eventually the core becomes un- 
stable and collapses; this is the start of the Class phase. At 
this stage the radius of the core is ~ 0.06 pc. As the subse- 
quent collapse proceeds, a freefall velocity field (v oc r _1//2 ) 
develops from the inside out, steadily replacing the more 
modest approximately uniform velocity field in the outer 
layers (the fourth and fifth rows on Fig. |l]) . The rate of ac- 
cretion onto the central protostar (sink) is quite modest, 
with a maximum value ~ 0.65 x 10 -5 Mq yr _1 , after which 
it decreases steadily into the Class I phase. The evolution of 
the sink mass and the accretion rate is shown on the first row 
of Fig. [?]. The durations of the prestellar and Class phases 
are given in Table 1, along with an estimate of the mean 
cruising velocity set up by the passage of the compression 
wave. 

This case is the one closest to the standard case of Foster 
& Chevalier (1993). In their simulation collapse was initiated 
by taking a marginally unstable (i.e. critical) Bonnor-Ebert 
sphere and increasing the density by 10%. As a result, con- 
traction was driven by an imbalance between self-gravity 
and pressure. Immediately following the formation of the 
central protostar the accretion rate was very high, and sub- 
sequently it decreased. In the <f> = 10 case here, contraction 
is driven by compression, and when the resulting compres- 
sion wave converges on the centre the accretion rate builds 
up to a maximum somewhat more slowly (than in Foster and 
Chevalier's simulation). In both cases the maximum accre- 
tion rate occurs during the Class phase, and is significantly 
larger than the standard M ~ a /G ~2x 1O" 6 M year -1 . 

Figure ^| shows the case <fr = 3 (i.e. quite slow compres- 
sion; P cxt doubles in 3 sound-crossing times) . The evolution 
is somewhat less quasistatic (than for (j> = 10, Fig. [l]), but 
basically it is very similar. The outer boundary is pushed 
inwards, and the core becomes steadily more centrally con- 
densed throughout the prestellar phase (the first three rows 
of Fig. ^); the outer layers start to move inwards at speeds 
in the range (0.05, 0.07) kms" 1 , whilst matter near the cen- 
tre is approximately stationary. Eventually the core becomes 
unstable, and the sink starts to accrete matter, marking the 
formation of the protostar and the start of the Class phase. 
At this stage the core radius is ~ 0.05 pc. A freefall velocity 
field then develops, as material accretes onto the protostar 
(fourth and fifth rows of Fig. H). The accretion rate reaches 
a maximum of ~ 1.08 x 10 M© yr _1 , and then declines 
monotonically (second row of Fig. r7j) . 

In Figure |^ we show results for = 1. Here the pressure 
doubles in one sound-crossing time, and so the evolution is 
more dynamic (than in the <fr = 10 and (j> — 3 cases de- 
scribed above). A small compression wave is driven into the 
core increasing the density and leaving in its wake a modest 
inward velocity field in the range (0.10, 0.15) kms -1 . The 
prestellar phase ends - and the Class phase begins - when 
this compression wave impinges on the centre. Up until this 
stage, the central density has hardly changed, since the inner 
parts have been unaware of the increased external pressure; 
the radius of the core has decreased to ~ 0.045 pc. During 
the Class phase, a freefall velocity field develops around 
the central protostar, but the outer parts of the envelope are 



Protostellar Collapse Induced by Compression 5 



me 0.297 (1 6 yr) 




0.02 0.04 0.06 0.08 
radius (pc) 

me 0.595 ( 1 6 yr) 




0.02 0.04 0.06 0.08 
radius (pc) 

me 0.832 ( 1 6 yr) 




0.02 0.04 

radius (pc) 



0.06 



me 1.202 ( 1 O 6 yr) 




0.02 0.04 
radius (pc) 



0.06 



me 1.258 (1 6 yr) 




0.02 0.04 
radius (pc) 



0.06 




0.00 0.02 0.04 0.06 0.08 
radius (pc) 




0.00 0.02 0.04 0.06 0.08 
radius (pc) 




0.00 



0.02 0.04 
radius (pc) 



0.06 



1 50 [ 



c 100 



50 



0.00 



0.02 0.04 
radius (pc) 



0.06 





50 


cm" 


40 


CP 






30 


b 






20 






c 


'0 


CD 




~o 






0.00 



0.02 0.04 
radius (pc) 



0.06 



9.o : 

2 



6 - 



9.8 i 



-3.0 



-19 
-19 
-20 



0L_ 
-3.0 



-2.5 -2.0 -1.5 

og(radius) (pc) 



-1 .0 




-3.0 -2.5 -2.0 -1.5 -1.0 
og(radius) (pc) 




-2.5 -2.0 -1.5 -1.0 

og(radius) (pc) 




-2.5 -2.0 -1.5 -1.0 
og(radius) (pc) 




-2.5 -2.0 -1.5 -1.0 
og(radius) (pc) 



Figure 1. The lefthand column shows velocity profiles on a linear-linear scale, the central column shows density profiles on a linear-linear 
scale, and the righthand column shows density profiles on a log- log scale, for the case <fi = 10.0 (very subsonic compression), at three times 
during the prestellar phase and two in the Class phase: first row, t = 0.30 Myrs; second row, t = 0.60 Myrs; third row, t = 0.83 Myrs; 
fourth row, t = 1.20 Myrs; fifth row, t = 1.26 Myrs. The dashed curve in the first row represents the initial state. We note that with 
^K,^^™^™ ~ ^ — in— 19 „„ m — 3 +~ ~ v,,,™!^*. A^^u.r ^ o ^ in— 4 „ m — 3 
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Figure 2. As for Fig. |l| but for the case <f> = 3.0 (subsonic compression): first row, t = 0.24 Myrs; second row, t = 0.48 Myrs; third row, 
t = 0.66 Myrs; fourth row, t = 0.80 Myrs; fifth row, ( = 0.87 Myrs. 
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Figure 3. As for Fig. |lj, but for the case </> = 1.0 (approximately sonic compression): first row, t = 0.24 Myrs; second row, t = 0.46 My] 
third row, t = 0.55 Myrs; fourth row, t = 0.59 Myrs; fifth row, t = 0.64 Myrs. 
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still moving inwards at approximately uniform sonic speed 
(v ~ 0.12 kms -1 ). The accretion rate is significantly higher 
than for the more quasistatic cases (larger </)), reaching a 
maximum of ~ 1.50 x 10 -5 MQyr -1 , and then decreasing 
into the Class I phase (third row of Fig. ^) . 

Figure ^ shows detailed results for the case <fi = 0.3 , i.e. 
quite rapid compression. In this case a strong compression 
wave is driven into the core, leaving a marginally sitteonic ve- 
locity field in its wake, i.e. v in the range (0.12, 0.14) kms -1 ). 
This leads to quite rapid accretion during the Class phase, 
reaching a maximum of ~ 1.65 x 10 -S M Q vr -1 , and then 
decreasing into the Class phase (fouth and fifth rows of 
Fig.fl). 

Figure fcd shows detailed results for the case <j> = 0.1 , 
i.e. very rapid compression. In this case a very strong com- 
pression wave is driven into the core, leaving a mildly 
supersonic velocity field in its wake, i.e. v in the range 
(0.25,0.30) kms" 1 ). At t = 0.24 Myr, the outer layers have 
been swept up into a tsunami. The density at the head of 
the tsunami is ~ 3 times higher than at the centre. This 
leads to very rapid accretion during the Class phase, reach- 
ing a maximum of ~ 2.60 x 1O _J M0 yr _1 , and then de- 
creasing into the Class I phase (fifth row of Fig. Q). In the 
Class phase, a freefall velocity field develops in the cen- 
tre, whilst the outer layers continue to cruise inwards at 
v ~ 0.30 kms -1 . 

If the external pressure only increases for a finite time, 
and then stays constant at a new higher value, the results 
are somewhat changed, particularly for the cases of rapid 
compression (small <j>). This is illustrated on Fig. ^, where 
we show detailed results for the cj> = 0.1 case when the pres- 
sure increase is halted as soon as the pressure has doubled. 
We refer to this case as l '<j> = 0.1 (finite)", to distinguish 
it from the standard case "(/> = 0.1 (indefinite)". In the case 
a (j> = 0.1 (finite)", the compression wave is much weaker and 
slower than for the case "0 = 0.1 (indefinite)", and so the 
inward velocity field it engenders is also weaker. In fact, the 
evolution for case "tj> — 0.1 (finite)" is more like that for the 
case "(/> = 0.3 (indefinite)" ; in particular, the compression 
wave takes a comparable time to converge onto the centre. 
For this reason, the mass and accretion rate for the case 
u (j> = 0.1 (finite)" are displayed with dashed lines on the 
same plot as the results for the case "0 = 0.3 (indefinite)" 
(i.e. the fourth row of Fig. Q). For slower compression rates, 
tf> <; 1 , halting the pressure increase when the pressure has 
doubled has little effect, because instability is triggered be- 
fore the pressure doubles, and the collapse dynamics of the 
core are therefore already well established. 

To test for convergence, we have repeated the case "0 = 
1 (indefinite)" with only A/lnt = 20, 000 particles in the core. 
The time-dependence of the protostellar mass and accretion 
rate are shown on the third row of Fig. (?| with dotted lines. 
We see that there is almost no change, except that with 
fewer particles the compression wave takes a little longer to 
converge on the centre, about 3% longer. We do not know 
the reason for this, but in any case it is a very small effect. 
Apart from this time delay, the detailed density and velocity 
profiles are essentially unchanged. 

Figure [7] shows the accretion rate onto the sink, as a 
function of time, for the various values of cf>. The star in- 
dicates where the Class phase gives way to the Class I 
phase. Table 1 shows the durations of the prestellar and 



Table 1. The durations of the prestellar and Class phases 
{Atprestellar an d Atci a sso)i an d the mean magnitude of the 
inward velocity field set up during the prestellar phase phases, 
^cruise : f° r different values of </> . 



<t> 


pre stellar 


A*ClassO 


^cruise 




(Myrs) 


(Myrs) 


(km S — 1 ) 


10.0 


1.2 


0.11 


0.025 


3.0 


0.8 


0.05 


0.06 


1.0 


0.6 


0.04 


0.12 


0.3 


0.4 


0.03 


0.17 


0.1 


0.3 


0.02 


0.30 



Class phases, for representative values of <j> . Since the sink 
particle is present from the outset, the moment of protostar 
formation is not precisely defined, and therefore these du- 
rations are not precisely defined. They simply illustrate the 
fact that the prestellar and Class phases are accelerated 
by rapid compression. In all cases, AtciassO / ' At preste iiar ~ 
0.05 to 0.10. 

We also tabulate values of demise , which is the mean in- 
ward velocity established in the outer parts of the core as a 
result of the passage of the inward propagating compression 
wave. These values are obtained from the detailed veloc- 
ity profiles by inspection, and are therefore only indicative 
of the general trend. The inward velocities observed in the 
outer envelopes of isolated prestellar cores in nearby star for- 
mation regions range from ~ 0.04kms~ 1 to ~ O.lOkms -1 
(Lee, Myers & Tafalla 1999). Our models generate compara- 
ble velocities for <j> in the range (1, 10), i.e. for situations in 
which the external pressure increases on a time-scale com- 
parable with or greater than the sound-crossing time. 

We speculate that collapses with lower values of <f>, i.e. 
more dynamical collapses, arise in more violently triggered 
star formation regions such as the p Ophiuchi and Perseus 
protoclusters (cf. Motte & Andre 2001). This would explain 
their having higher absolute densities and smaller extents at 
the onset of the Class phase. 



4 DISCUSSION 

The obvious inference to be made from these simulations 
is the unsurprising one that, the more rapid and the more 
prolonged the increase in external pressure, the more rapid 
is the evolution, the stronger is the compression wave, the 
higher is the peak accretion rate, and the shorter are the 
prestellar and Class phases. However, other important in- 
ferences are to be found in the detail. In particular, simula- 
tions with <j> <; 1 (i.e. pressure increasing on a sound-crossing 
time or slower) match essentially all the observational con- 
straints described in Section 2. 

First, the density profiles in the prestellar phase are 
relatively flat in the inner regions, r/ = —d£n[p]/d£n[r] ^ 1 , 
and steepen towards 77 ~ 2 in the outer parts. The effect of 
the compression wave will be hard to detect in projection. It 
will simply produce an apparent extension of the relatively 
flat inner part of a prestellar core. By the Class and Class I 
phases, the density profiles should steepen, viz. to r\ ~ 3/2 in 
the inner regions, and to r\ ~ 2 — 3 in the outermost regions 
(see log- log plots of the density at late times in Figs, [j] to ^[ 

Second, the velocity fields are very reminiscent of that 
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Figure 4. As for Fig. |lj, but for the case </> = 0.3 (mildly supersonic compression): first row, t = 0.12Myrs; second row, t = 0.24Myi 
third row, t = 0.36 Myrs; fourth row, t = 0.43 Myrs; fifth row, t = 0.46 Myrs. 
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Figure 5. As for Fig. 
row, t = 0.32 Myrs; fift 



lj, but for the case c/> = 0.1: first row, t = 0.12 Myrs; second row, t = 0.24 Myrs; third row, t = 0.29 Myrs; fourth 
l row, t = 0.34 Myrs. 
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Figure 6. As for Fig. |lj, but for the case <f> = 0.1 (finite), where instead of letting P e xt increase indefinitely, the increase is halted once 
P ext has doubled: first row, t = 0.12 Myrs; second row, t = 0.24 Myrs; third row, t = 0.33 Myrs; fourth row, t = 0.42 Myrs; fifth row, 
t = 0.46 Myrs. 
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Figure 7. Accretion rates as a function of time. First row, <^> = 10 ; second row, <j> = 3 ; third row, <j> = 1 ; fourth row, <j> = 0.3 ; fifth row, 
(f> = 0.1 . The full curves correspond to the standard simulations, i.e. those in which P e xt is allowed to increase indefinitely. The dotted 
curves on the third row correspond to the "<j> = 1 (indefinite)" case, but performed with fewer particles (i.e. A/i n t = 20, 000 instead of the 
standard M ln t = 50, 000) in the core. The dashed curves on the fourth row correspond to the case <^> = 0.1 (very fast compression), but 
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inferred for L1544 by Tafalla et al. (1998) and Williams et 
al. (1999). In L1544 the velocity field in the outer parts of 
the core appears to be approximately uniform, with magni- 
tude ~ 0.08 km s -1 . Other cores showing similar infall sig- 
natures give inward velocities in the range ~ 0.04 km s -1 to 
~ 0.10 kms" 1 (Lee, Myers & Tafalla 1999). If these velocities 
are due to an inward-driven compression wave of the type 
we have simulated, it suggests that the external pressure 
increases on a time-scale comparable with, but somewhat 
greater than, a sound-crossing time. 

The generation of an approximately uniform inward ra- 
dial velocity field appears to be a generic feature of com- 
pression waves. All of the two-dimensional infinity of simi- 
larity solutions found by Whitworth & Summers (1985) in- 
volve compression waves which are driven from the outside 
inwards, and which leave behind them a uniform inward 
velocity field (v ra d — constant) and a p oc r~ 2 density 
field. The convergence of this compression wave on the cen- 
tre at t = signals the formation of a central stellar ob- 
ject. This object then grows at a constant accretion rate. 
At the same time an expansion wave is reflected outwards 
and leaves in its wake an approximately freefall velocity 
field, -u rad oc r~ 1/2 , and p oc r' 3/2 . The Shu (1977) solu- 
tion is simply the limiting case of a zero-amplitude com- 
pression wave leaving a zero-amplitude inward velocity field 
in its wake; hence, uniquely, it starts from a static singular 
isothermal sphere. The Larson-Penston solution involves a 
fairly strong compression wave propagating into a uniform- 
density, homologously contracting cloud. If radial velocity 
fields like that inferred for L1544 turn out to be common, 
one possible explanation is that core collapse and protostar 
formation generally follow from external compression. 

Third, the accretion rate is very high immediately after 
the formation of the central protostar, i.e. in the Class 
phase, but later it decreases. This conforms to the pattern of 
accretion inferred from the decline in outflow power between 
the Class and the Class I phase (Bontemps et al. 1996, 
Henriksen et al. 1997), and is consistent with the relative 
numbers of Class and Class I sources. 

A rapid increase of external pressure, of the sort we 
have invoked in these simulations, could arise in a number 
of situations. In the first instance, we might presuppose the 
existence of an approximately stable gravitationally bound 
core. If this core were then overrun by a strong shock which 
propagated supersonically through the medium surrounding 
the core, the external pressure acting on the core bound- 
ary would be increased. Alternatively, if the medium sur- 
rounding the core were suddenly irradiated more intensely, 
because a luminous star switched on nearby, or came out 
from behind a dust cloud, the surrounding medium would be 
heated - and possibly also ionized - giving a rapid increase in 
the external pressure. A contraction of the molecular cloud 
in which the core is embedded, as in the scenario of Palla & 
Stahler (2000), would also lead to an increase (presumably 
slow) in the external pressure bounding the dense core. A 
similar type of collapse would ensue if the radiation heat- 
ing the core were suddenly shut off (but not that heating 
the external medium). The core would become thermally 
unstable, and the external medium would drive in a com- 
pression wave. It is also possible that the same convergent 
flows which create cores in the first place, may sometimes 
then continue to compress these cores until they collapse, as 



in the scheme recently described by Hartmann, Ballesteros- 
Paredes & Bergin (2001). 

However, it is very unlikely that the increase in the ex- 
ternal pressure would actually be isotropic - as we have 
assumed, for simplicity, in the simulations presented here. 
Therefore we should not expect compression waves to con- 
verge on the centre simultaneously from all directions. In- 
stead, the compression waves arriving from different direc- 
tions are likely to be out of phase, so that they will interfere 
in complex ways. This could have the effect of creating mul- 
tiple protostellar seeds in the interior of a core, in accordance 
with the observation that most stars are born in binary or 
higher multiple systems (Mathieu 1994). 



5 CONCLUSIONS 

We conclude that the model of protostellar collapse devel- 
oped here merits further investigation. It appears to repro- 
duce the gross features of the density and velocity fields 
observed in prestellar cores and protostars, and the relative 
ages of the Class and Class I phases. We now need to ex- 
plore the consequences of introducing an energy equation, 
the effect of initial rotation, and the influence of perturba- 
tions as a means of inducing the formation of binaries and 
higher multiples. 
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